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Abstract. For a sample of 51 blazars with extensive optical polariza- 
tion data, we used circular statistics to calculate the scatter among the 
polarization position angles for each object. We found that this scatter is 
correlated with the radio core dominance. We compared this relationship 
with the predictions of a simple transverse shock model. The result sug- 
gests that blazar jets are likely to cover a wide range of speeds consistent 
with those derived from observations of superluminal motion. 



1. Introduction 

Jet are very important for AGN unification. All radio loud objects have powerful 
jets. In the orientation Unification Scheme, blazars, with their extreme flux and 
polarization variability, have the jets beamed towards us. However viewing 
angle unification is not the complete story. Do all the jets have similar physical 
properties on average, so that orientation is able to explain most of the variations 
of observed properties among blazars? 

Optical polarization of blazars provides magnetic field information very 
close to the jet formation region and is therefore a very important tool for 
probing jet physics. We have investigated a sample of 51 blazars with many 
optical polarization measurements made over the past 20-30 years. In paper I, 
we have investigated the relationship between the optical polarization position 
angle and the VLBI radio structure. The results are consistent with a shock 
model. In this paper, we focus on the variability of optical polarization position 
angles. Assuming a simple shock model, we link the angle variability with jet 
orientation and speed to explore the jet physical parameters among blazars. 



2. Data and Analysis 
2.1. Optical Polarization 

Since a blazar's optical polarization swings over a wide range of angles, there 
exists a 180° ambiguity. Linear statistics based on normal distributions is not 
strictly correct. Therefore we use circular statistics - a technique little used in 
astronomy but commonly used for example, to analyze directions of returning 
homing pigeons, or wind directions. Thus, from the polarization position angles 
9i, we calculate the circular standard deviation c(9) about the mean angle, to 
represent their scatter (Fisher 1993). Polarization data are axial, i.e. only cover 
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the range 0° to 180°. However, the 29 i are truly circular data. Therefore, we 
first calculate c(29), then divide by 2 to get c(9): 



c{9) 



ln ( V(EF=i cos 29 t y + (E-'=i sin 2^ 



where n is the total number of measurements. If the sample has a large intrinsic 
scatter, the determination of the real scatter depends on the accurate shape of 
the wrapped distribution between —tt/2 and tt/2. If the sampling of data is 
not enough to reveal the detailed shape of the distribution, the c(9) calculation 
algorithm gives a maximum c(9) ~ 66°. 

2.2. Radio Core Dominance 

In the unification scheme of powerful radio sources, the core-dominated sources 
are the Doppler-boosted jet, viewed end-on. At larger angles the core becomes 
much weaker and collimated jets link the well resolved double lobes. Doppler 
boosting in core-dominant sources is supported by their very high brightness 
temperatures and by superluminal motions. Statistical analyses of radio source 
samples support the unified scheme (M. Chiaberge, this conference). The ratio of 
core-to-lobe luminosity, or the core dominance, R, is therefore a good indicator 
of viewing angle. However, in blazar samples, viewing angle is not the only 
parameter affecting Doppler boosting. From work of other authors (Orr and 
Browne 1982), we know that R depends on three physical parameters and can 
be expressed as 

logR = log(R T x (1 - x cosO)) -3 ) 

where <f) is the viewing angle of the jet, (3 is the bulk speed of the emitting 
electrons in the jet and Rt is the tangential value of R measured for an average 
edge-on jet in radio galaxies. We used Rt = 0.0024 at 5GHz rest frequency (e.g. 
Figure 1 of Hoekstra, Barthel & Hes 1997). While R is of course sensitive to <p, 
for smaller <f> (e.g. < 20°), it is increasingly sensitive to (3. 

For each object we calculated R at 5 GHz rest frequency, using core and 
lobe flux densities from literature. In Section 4, we will calculate logR from 
models using the formula above and compare the models with observed values. 



3. Correlation 

Figure 1 (a) shows that there is a correlation between c{9) and log R. The corre- 
lation is significant at the 99.9% level. However, we note that the correlation is 
not linear. For the logR<2 region, which includes most blazars, c(9) increases 
with logR, but for logR>2, there is a lot of scatter. No object has c(9) larger 
than 66°, as expected from the algorithms for large scatter and small sample 
size. 
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Figure 1. (a) - Correlation between log R and the scatter of optical 
polarization position angle. Triangles are BL Lacs and squares are 
QSOs. The two thin curves represent (3 = 0.93, 0.99 models, (b) - 
Model simulations. Each represents represents a model with fixed (5 
and variable <f>. The (5 for each model is noted on top of each curve. The 
markers along each curve denote the corresponding (j) in the observer's 
frame. 



4. Model 

The model we use makes the following assumptions: 

1. Optical emission comes from a transverse shock propagating along the jet. 

2. The shock compresses the magnetic field so much that the magnetic com- 
ponent along the jet direction is negligible. 

3. In the plane of the shock, there is a dominant magnetic field direction that 
changes randomly with time. This magnetic field direction corresponds 
to the emission region whose radiation happens to be boosted towards us. 
Different ejection directions observed for emerging blobs on VLBI scales 
and complex VLBI structure suggests that these dominant regions change 
rapidly. 

4. The shock is optically thin. 

When the shock is viewed face on, the magnetic field direction projected 
on the sky will vary randomly with time. As the viewing angle is increased, 
the projected magnetic field direction will lie increasingly perpendicular to the 
projected jet direction. 

A simplification is afforded by the Lorentz invariance of the ratios of Stokes 
parameters, so our calculations of projected magnetic field direction in the co- 
moving frame (viewing angle (j)') can be referred to the observer's frame using 
the relation (Bjornsson, C.-I. 1982): 
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1 — pcos((p) v 

For a given set of (3 and 0, we generate 100 randomly oriented magnetic 
field directions to simulate the random changes over time inside the shock plane 
and project them onto the sky in the co-moving frame. The circular standard 
deviation of those angles can then be compared with real observations, using 
the relation between R and 4> given earlier. Rt is fixed by observations, so the 
only free parameter in this model is (3. 

The results are shown in Figure 1 (b). Each curve represents a model with 
fixed (3, with <p decreasing from bottom left to top right. The values of 4> are 
marked along the curve. From left to right, (3 = 0.93,0.95,0.97,0.99. The 
upper limit on c{9) is the result of the computational limit. The lower limits are 
caused by the fact that we only have 100 simulated vectors to do statistics and 
any one vector that happens to be parallel to the jet direction contributes a lot 
to the scatter. In principle, if we had an increasing number of vectors, the lower 
limit for c(9) would approach 0. On our model, the c(0) vs. logR correlation 
arises because (3 ~ 0.95 for most blazars, and so the viewing angle dependence 
dominates for logR< 2. For logR> 2, the scatter is the result of a tail to higher 
jet speeds. 

Our simple model with only one free parameter, (3, accounts for the cor- 
relation we find between the scatter in optical polarization angles, c(9), and 
the core-dominance, log R. Figure 1 shows that the average jet speed differs 
significantly amongst the powerful blazars, with a range of values comparable 
with those derived from superluminal motion and from the statistics of core- 
dominance in radio-source surveys (Urry & Padovani 1995). In Unified Schemes 
it is important to take into account that log R is not simply a measure of orien- 
tation. Jet speed introduces considerable scatter, especially at log R > 2. 
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